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ABSTRACT 

The X-ray burster GS 1826-24 shows extremely regular Type I X-ray bursts whose energetics and recur- 
rence times agree well with thermonuclear ignition models. We present calculations of sequences of burst 
lightcurves using multizone models which follow the nucleosynthesis (ap and rp-processes) with an extensive 
nuclear reaction network. The theoretical and observed burst lightcurves show remarkable agreement. The 
models naturally explain the slow rise (duration w 5 s) and long tails 100 s) of these bursts, as well as 
their dependence on mass accretion rate. This comparison provides further evidence for solar metallicity in the 
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accreted material in this source, and gives a distance to the source of (6.07 ±0.18) kpc where is the 

burst emission anisotropy factor The main difference is that the observed lightcurves do not show the distinct 
two-stage rise of the models. This may reflect the time for burning to spread over the stellar surface, or may 
indicate that our treatment of heat transport or nuclear physics needs to be revised. The trends in burst proper- 
ties with accretion rate are well-reproduced by our spherically symmetric models which include chemical and 
thermal inertia from the ashes of previous bursts. Changes in the covering fraction of the accreted fuel are not 
required. 

Subject headings: accretion, accretion disks — stars: individual (Ginga 1826-238, GS 1826-24) — stars: neu- 
tron — X-rays: bursts 



1. EVITRODUCTION 

The basic physics of Type 1 X-ray bursts as thin shell 
flashes on the surfaces of accreting neutron stars was under- 
stood many years ago (e.g. lFujimoto et al]ll981tlLewin et alJ 
Il995h . Nonetheless, detailed comparisons of observations and 
theory are often less th an successful (Fujimoto et al. 1987; 
Ivan Paradij set al.lll988l) . The burster GS 18 26-24 (aka Ginga 
1826-24) is an important exception. Ubertini et al.l (Il999h 
dubbed it the "clocked burster" because of its extremely reg- 
ular bursting behavior. They found a burst recurrence time 
close to 6 hours with a dispersion of approximately 6 min- 
utes. iBildsteni (2000) dubbed it the "textbook burster" be- 
cause of the good agreement with theory. He noted that 
the burst energetics and recurrence times were as expected 
for the inferred accretion rate M « 10"' Mp) yr"\ and pro- 



posed that the long burst tails, last ing ; 



yr 

100 s, were po wered 



by rp-process hydrogen b urning (I Wallace & Wooslevlll98ll : 
iHanawa & Fujimotolll984l) . The rp-process involves a series 
of proton captures and beta-decays on heavy nuclei close to 
the proton drip line, for which there are considerable uncer- 
tainties in beta decay rates, reaction rates, and nuclear masses. 
Type 1 X-ray bursts offer an important test of our understand- 

in g of this process . 

iGallowav etall (|2004 (hereafter G04) studied 24 bursts 
observed from GS 1826-24 by the Rossi X-ray Timing Ex- 
plorer (RXTE) between 1997 November and 2002 July, and 
carried out a detailed comparison of the observed recurrence 
times and energetics with theory. During this period, the 
accretion rate (assumed to be proportional to the persistent 
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X-ray luminosity) increased by « 50%, while the burst re- 
currence time decreased from Af « 6 hours to « 4 hours. 
At the observed accretion rate M « 10"' Mq yr"', hydrogen 
bur ns stably between bursts by the beta-limited hot CNO cy- 
cle (iHoyle & Fowled ri965i) . heating the accumulating layer 
of hydrogen and heUum at a rate that depends only on the 
mass fraction of CNO elements. The resulting ignition col- 
umn de2tii_of_~^_2<_20^cnr^ is nearly independent of M 
(e.g. lCumming & Bildstenll2000h . in agreement with the ob- 
served scaling which is close to Af c>c 1 /M. 

G04 used the variation of burst energetics with accretion 
rate to constrain the composition of the accreted material. The 
energetics are measured by the ratio a of integrated persistent 
flux between bursts to the burst fluence, equivalent to the ra- 
tio of the gravitational energy release from accretion to the 
nuclear energy release in bursts. The observed value a « 40 
is close to the value expected for the 4-5 MeV per nucleon en- 
ergy release in hydrogen burning (Bildsten 2000). G04 found 
that a decreased from « 44 in 1997 to « 40 in 2002. This is 
consistent with the change in composition of the fuel layer ex- 
pected from hot CNO burning between bursts, if the accreted 
material has solar metallicity. For a shorter recurrence time, 
less hydrogen burns between bursts, resulting in a larger burst 
energy and smaller a. In addition, G04 found that whereas 
the ignition models predicted that the ignition column depth 
should increase slightly with M, the observed trend was the 
opposite. 

The study of G04 did not include models of the burst 
lightcurves. In this paper, we cor npare observations o f 
GS 1826-24 with multizone models dWooslev et all |2004|) . 
which include a large nuclear reaction network to follow the 
rp-process in detail. We show that these models naturally ex- 
plain the slow rise (duration w 5 s) and long tails (w 100 s) of 
these bursts. We show that our time-dependent models repro- 
duce the scalings of the recurrence time with accretion rate 
for a metallicity near solar. 
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Fig. 1. — Comparison of observed and calculated lightcurves. The his- 
togram shows the average Hghtcurve from the bursts observed during the year 
2000 when the recurrence time was fa 4 hours (G04. Figure 2). The error bars 
are the Icr variations from burst to burst. The solid and dashed curves are the 
average burst profiles from models A3 (Z = 0.02) and B3 (Z = 0.001), which 
have A; = 3.9 and 4.0 hours respectively. The inset magnifies the rise and the 
initial part of the decay. The grey bands indicate the 1 u variation of the burst 
profiles about the average. 

2. COMPARISON OF OBSERVATIONS WITH TIME-DEPENDENT 
BURST MODELS 

We calcula te sequences of bursts as described by 
IWooslev et ail ([2004) (hereafter W04). The KEPLER code 
dWeaver et al. [l978) is used to follow the evolution of the 
outer layers of the neutron star through a sequence of succes- 
sive bursts. A large adaptive nuclear reaction network is used 
to follow the nucleosynthesis at each depth, and we include 
convection when needed using time-dependent mixing length 
theory. We use the same input nuclear physics, stellar opaci- 
ties, and neutron star parameters as W04, but consider a wider 
range of accretion rates. The Newtonian calculations are cor- 
rected for general relativity as described in §4.4 of W04 for a 
neutron star mass of 1 .4 Mq . The corresponding stellar radius 
is 7? = 11 .2 km, and gravitational redshift is z = 0.26. 

A summary of the results is given in Table [T] For each 
sequence of bursts, we list the rest mass accretion rate M, 
recurrence time Af, burst energy fiburst, gravitational mass 
accumulated between bursts AM = MAf/(l+z), and a = 
FxAf/iiburst = AMc^z/iiburst, all as seen by an observer at in- 
finity. The quantities given are averaged over all bursts except 
the first burst in each sequence, which is typically more en- 
ergetic than the subsequent bursts (W04). We give the stan- 
dard deviation of each quantity in parentheses, to show the 
burst to burst variations. Models A4 and B3, which have 
M = 1.75 X 10"^ Mq yr~', have the same parameters as mod- 
els ZM and zM of W04. There are slight differences at the 
level of w 3% between the burst properties in models A4 and 
B3 as compared to models ZM and zM of W04, because of 
refinements of the KEPLER code that were made following 
publication of the W04 paper 

We compare these simulations to bursts observed by RXTE 
between 1997 November and 2002 July. We analyze the data 
as described in G04, with the following exceptions: 1) the 
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Fig. 2. — Burst pai'ameters for the Z = 0.02 models A1-A4 compared with 
individual observed bursts. The properties of the observed bursts are plotted 
as diamonds. The properties of individual model bursts are plotted as crosses, 
with average values for each model plotted as squares joined by solid lines. 

spectral fitting was performed using LHEASOFT version 5.3, 
released 2003 November 17, for which the effective area of 
the proportional counter instrument (and hence the source 
flux) was reduced by approximately 15% compared to earlier 
versions; 2) improved calculation of the burst fluence which 
was better able to handle gaps in the data, which increased 
the estimated fluence in some cases by at most 5%. These 
changes also had the effect of reducing the absolute a-values, 
although the trend with persistent flux was unchanged. 

In Figure [T] we compare the mean Hghtcurve for bursts 
observed during 2000 (G04 Figure 2) with the mean burst 
lightcurves from models A3 and B3. These models are chosen 
because they have similar recurrence times to the observed re- 
currence time of 4. 1 hours. We calculate the mean lightcurves 
by aligning bursts in each sequence by their peak luminosities. 
The error bars in Figure[T]show the la burst to burst variation 
about the mean observed Hghtcurve. The shaded region shows 
the same variation for the theoretical lightcurves. 

For this comparison, we choose the distance to the source 
(within the allowed range 4 < c/ < 8 kpc; G04) so that the 
peak luminosity of the observed bursts agrees with the peak 
luminosity of bursts from model A3. The relation between 
the peak burst luminosity Lpeak and the observed peak flux is 
47rii^^iFpeak = ipeak, wHcre (^h is a factor that accounts for pos- 
sible anisotropy in the burst emission (e.g. Fujimoto 1988). 
The average observed peak flux in the 2000 epoch is (2.93 ± 
0. 15) X 10"*^ erg cm"^ s"' , and the average peak luminosity of 
bursts in model A3 is Lpeak = (1-29 ±0.04) x 10^^** erg s"', giv- 
ing a distance d = (6.07 ±0.18) kpc f^^^^'. 

Once the bursts have been normalized in this way, the 
agreement between the observed and theoretical lightcurves 
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TABLE 1 
Average Burst Properties" 



Model Number Z M At Ei,„st a AM 

of bursts (10-' M0 yr-l) (h) (lO^' ergs) (10^' g) 



Al 
A2 
A3 
A4 
Bl 
B2 
B3 



19 
18 
30 
13 
12 
17 
15 



0.02 
0.02 
0.02 
0.02 
0.001 
0.001 
0.001 



1.17 
1.43 
1.58 
1.75 
1.17 
1.43 
1.75 



5.4 (0.1) 
4.3 (0.1) 
3.85 (0.06) 
3.48 (0.06) 
12.8 (0.6) 
6.04(0.41) 
3.98 (0.28) 



4.67 (0.20) 
4.67(0.11) 

4.73 (0.07) 
4.84 (0.06) 

13.3 (0.7) 

7.74 (0.49) 
6.26 (0.32) 



57.4 (2.8) 
55.6(1.2) 
55.0(0.9) 
53.6 (0.8) 
47.8 (0.4) 
47.4(1.0) 
47.3 (2.4) 



1.14(0.03) 
1.11 (0.03) 
1.10(0.02) 
1.11 (0.02) 
2.71 (0.14) 
1.57 (0.11) 
1.27 (0.09) 



' See text for definitions of quantities. 
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Fig. 3. — Change in the lightcurve as a function of accretion rate. Averaged 
lightcurves from models Al (dotted curve), A2 (solid curve), and A3 (dashed 
curve). The grey shaded areas indicate the standard deviation of the burst to 
burst variation ai'ound the average lightcurve. As accretion rate increases, the 
burst becomes less luminous, with a longer rise time and slower decay. 

for Z = 0.02 (model A3) is remai-kable. Model A3 fits the 
observed decay exceedingly well out to about 40s, falling a 
little below the observed flux between 40 and 120s. The burst 
to burst variations in the models are also of comparable mag- 
nitude to the burst to burst variations in the data. The most 
significant difference is that the theoretical model shows a 
distinct two-stage rise which is not apparent in the observed 
lightcurve (see the lower panel of Fig. 1). 

Model B3, which has a low metallicity, does not reproduce 
the observed lightcurve. Given the uncertainty in the distance 
to the source, the normalization may be adjusted to bring the 
observed and predicted peak luminosities into agreement, but 
the shape of the decay provides an additional constraint. In 
model B3, the lower metallicity leads to very little hydrogen 
burning between bursts, giving less helium and more hydro- 
gen at ignition than in the solar metallicity model. The result 
is a burst with a much longer tail, inconsistent with the ob- 
served profile for any distance. 

In Figure |2l we compare the burst recurrence times, ener- 
gies, a values, and ignition masses as a function of accretion 
rate for the observed bursts and the solar metallicity models 
A1-A4. The observed burst fluence Eb is related to the burst 
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Fig. 4. — Fitted exponential decay timescales for observed bursts (open 
symbols) and theoretical lightcurves from models Al, A2, A3 and A4 (solid 
symbols). For the models, we choose a representative sample of 5 bursts 
(burst numbers 7 to 11 from each run), and displace them slightly on the 
horizontal axis for clarity. Note that some of the theoretical lightcurves show 
a three stage decay. The uncertainties in the fitted decay times are comparable 
to the scatter in the points. We adopt the same relation between Fx and M as 
in Figure |2] 

energy fiburst by Ei, = A-nd^ £^bEhuKt, and we adopt the same 
value of (_hd^ that we determined earlier by comparing the 
burst peak luminosity with the observed peak flux. We convert 
the bolometric flux Fx to accretion rate M using the relation 
4-TTd^^pFx = Mc^z/(l +z), where is a factor that accounts 
for possible anisotropy in the persistent emission (e.g. Fu- 
jimoto 1988). For the same value of distance, we find that 
£,p/£,b = 1-55 is required for the recurrence times to match in 
Figure |2] The theoretical alpha values have been divided by 
this same factor to give the correct prediction for alpha includ- 
ing anisotropic emission. 

The agreement between the models and observations in Fig- 
ure|2]is excellent. The main difference is that the burst energy 
is overpredicted by approximately 5%, and the a value un- 
derpredicted by the same factor The low metallicity models 
B1-B3 are not consistent with the observed trends. For ex- 
ample, the burst energies show a factor of 2 variation in these 
models over the observed accretion rate range (Table 1), much 
larger than the « 10% variations seen in the data. 
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Figure [3] shows how the burst Hghtcurves change with ac- 
cretion rate (models Al to A3). G04 pointed out that the 
burst Hghtcurves changed with recurrence time. They fit- 
ted the decay of the bursts to an exponential profile with a 
break in the exponential decay time. The theoretical curves 
show a similar two timescale decay (this is particularly ap- 
parent in the lower panel of Fig. 1). G04 noted a significant 
increase in the first exponential decay timescale as Af de- 
creased, from 14.7 ± 0.7s when Af w 6 hours to 17.5 ± 1.1s 
when Af ?» 4 hours. The second exponential decay timescale 
43 ± Is showed no significant change. G04 speculated that 
this change was a result of the composition of the fuel, since 
more hydrogen is depleted between bursts for larger Af . 

Figure [3] shows that the model Hghtcurves show a similar 
change with accretion rate to the observed Hghtcurves, pre- 
sumably relating to the changing hydrogen/helium ratio at 
ignition. This is quantified in Figure |4] which shows fitted 
exponential decay times for the observed and model bursts. 
Overall, the agreement is good, although the decay times 
of the model bursts are generally slightly shorter than ob- 
served (consistent with the difference in the burst tails in Fig- 
ure 1), and they show much larger burst to burst variations 
than observed. Also, whereas the observed bursts are well fit 
by a two-timescale decay, some of the model bursts show a 
more complex decay profile, which we have fitted with three 
timescales. 

3. DISCUSSION 

We have presented a first comparison between multizone 
burst models and the observed bursts from GS 1826-24. The 
good agreement further confirms GS 1826-24 as a "textbook" 
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burster. For a distance of (6.07 ±0.18) kpc (where ij, 
is the burst emission anisotropy factor), our solar metallic- 
ity models closely match the observed Hghtcurves. Both the 
long rise and decay times arise naturally from rp-process hy- 
drogen burning. From the very regular nature of the burst- 
ing, G04 argued that the accreted material covers the entire 
surface of the neutron star; this is supported by the excellent 
agreement of our spherically-symmetric models. Solar metal- 
licity models agree best with the observed Hghtcurves. Low 
metallicity models produce too little helium by hot CNO burn- 
ing prior to ignition, leading to a lower peak luminosity and 
a longer rp-process tail. This agrees with the conclusions of 
G04, who argued for solar metallicity based on the burst ener- 
getics. The estimate of the distance is based on a comparison 
of mean burst lightcurve at a single epoch, with the lightcurve 
predicted by the model for parameters giving the same recur- 
rence time. It is possible that lightcurve comparisons at dif- 
ferent epochs (i.e. values of the recurrence time) may result 
in different values of the distance and/or anisotropy parame- 
ters. Such comparisons will provide an additional consistency 
check for the distance estimation, which we will undertake in 
a future paper. 

Our models reproduce the observed variation in burst recur- 



rence times, energies, and a values with accretion rate for a 
ratio of anisotropy factors for the persistent and burst emis- 
sion of ^o/^fe = 1.55. This is within the range discussed by 
Fujimoto (1988) and similar to that inferred for other sources 
(e.g. Sztaino et al. 1987). G04 found that the ignition models 
of C umm in g & Bildsten (2000) reproduced the change in a, 
but not the scaHng of ignition column with M. In these mod- 
els, there is less time for helium production between bursts at 
higher M, delaying the ignition, and leading to increasing ig- 
nition mass. However, our time-dependent models show that 
the chemical and thermal inertia associated with the ashes 
from previous bursts (Taam 1980) outweighs the lower he- 
lium abundance, heating the layer, and leading to a smaller 
ignition mass as M increases, in agreement with observations. 
A change in the covering fraction of accreted fuel, as specu- 
lated by G04, is not required. 

Despite the good agreement, there are several differences 
between the observations and the models which need to be in- 
vestigated in future work. These comparisons promise to con- 
strain the nuclear physics of the rp-process, the thermal prop- 
erties of the burning layers, spreading of the nuclear burn- 
ing, and neutron star parameters. The model burst Hghtcurves 
show a distinct two-component rise which is not present in 
the data, as well as a slightly shorter burst tail. Another dif- 
ference in burst shape is that, unlike the observations, some 
of the model Hghtcurves show a three-stage rather than two- 
stage decay. These differences may relate to our treatment 
of heat transport (especially convection) or nuclear physics. 
W04 confirmed previous findings that b urst tails are sensi- 
tive to nuclear flo ws above the iron group dSchatz et alj2001al : 
iKoike et aLlll999 ). and also pointed out that the rise times are 
sensitive to nuclear decays below the iron group. An alterna- 
tive for the rise is that a finite propagation time for the burn- 
ing around the star of 1 s (e.g. Frvxell & Wooslev 198^ 
ISpitkovskv et al.li2002 ) might act to "wash out" the kink. 

The neutron star mass and radius also change the predicted 
burst properties. In this paper, we have considered a neutron 
star with redshift factor z = 0.26, corresponding to M = 1 .4 Mq 
and R= 1 1 .2 km. A smaller radius of 10.6 km would reduce 
the predicted burst energies by 5%, and would increase the 
redshift factor and therefore a, bringing both of these quanti- 
ties into better agreement with observations (Fig. 2). We will 
investigate the constraints on the neutron star mass and radius 
in detail in future work. 
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